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We report a scattering matrix theory for dynamic 
and nonlinear transport in coherent mesoscopic con- 
ductors. In general this theory allows predictions of 
low frequency linear dynamic conductance, as well 
as weakly nonlinear DC conductance. It satisfies the 
conditions of gauge invariance and electric current 
conservation, and can be put into a form suitable for 
numerical computation. Using this theory we exam- 
ine the third order weakly nonlinear DC conductance 
of a tunneling diode. 
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I. INTRODUCTION 

Quantum transport under a time-dependent field in co- 
herent jngsoscopic systems is the subject of many recent 
studiesErQ. Another problem of interest is the nonlinear 
conductance of such a_ system, whether under a time- 
dependent field or notfi. A difficult theoretical issue is 
the prediction, for a general mesoscopic conductor, the 
transport coefficients as a function of the AC field fre- 
quency and the bias voltage. Once these parameters are 
known, one can predict useful information such as the 
nonlinear current-voltage characteristics in the DC case, 
the emittance in the linear frequency, linear voltage AC 
case, and further nonlinear dynamic conductance. In- 
deed, it is now possible to experimentally measure the 
nonlinear AC transport properties such as the second 
harmonic generation, as have been demonstrated by sev- 
eral laboratoriesErtll. 

When a conductor is subjected to time-varying exter- 
nal fields such as an AC bias voltage, the total elec- 
tric current flowing through the conductor consists of 
the usual particle current plus the displacement current. 
The presence of the latter is crucial such that the total 
electric current is conserved. Hence for a theory to deal 
with AC transport, in principle one should include the 
displacement current into the consideration. Because a 
displacement current originates from induction, and the 
necessary condition for electric induction is the electron- 
electron (e-e) interaction, one thus concludes that an im- 
portant ingredient for AC transport theory should be 
the consideration of e-e interactions. These issues have 
been emphasized by Biittiker and co-workerslia On the 
other hand, for DC transport under nonlinear conditions, 
a necessary requirement is the gauge invarianceu: the 
physics should not change when electrostatic potential 



everywhere is changed by the same constant amount. 
Gauge invariance puts severe conditions on the form of 
the nonlinear transport coefficients. From these physical 
arguments, it is clear that AC as well DC nonlinear trans- 
port contains ingredients which were not needed when 
dealing with the familiar DC linear transport!! 2 ! 

The problems of current conservation and gauge invari- 
ance have been recognized in the literature. For conduc- 
tors which maintain quantum coherence, Biittiker and 
his co-workers have developed an approacho based on 
the single electron scattering matrix theory to deal with 
the linear AC dynamic conductance as well as the sec- 
ond order nonlinear conductance coefficients. The origi- 
nal scattering matrix theory were invented to investigate 
DC linear transport coefficients—as is represented by the 
Landauer-Biittiker formulationtj. Such a theory calcu- 
lates particle current from the scattering matrix, thus a 
direct application to AC situation would violate current 
conservation! 2 !!! 2 ]. To solve this problem, the scattering 
matrix theory for AC transport consists of two stepsEII!3. 
First, it calculates the particle current and finds that this 
current is not conserved. Second, it considers the e-e in- 
teraction which alters the scattering potential landscape, 
and this effect generates an internal response which can- 
cels exactly the non-conserved part of the particle cur- 
rent thereby restores the current conservation. For the 
DC second order nonlinear conductance coefficient, sim- 
ilar considerations led to the desired gauge invariance. 

In a recent work, the authors have developed a micro- 
scopic and general theoretical formalism for electric re- 
sponse which is appropriate fopboth DC and AC weakly 
nonlinear quantum transporto. That formalism was 
based on the response theory and it formalized the con- 
nection of the response theory to the scattering matrix 
theory at weakly nonlinear level. One of the useful con- 
ceptual advances of the general formalisrrJlil was the in- 
troduction of a frequency dependent characteristic poten- 
tial at the nonlinear level. The characteristic potential 
describes the changes of scattering potential landscape of 
a mesoscopic conductor when the electrochemical poten- 
tial of an electron reservoir is perturbed externallyn. It 
is the nonlinear order characteristic potential which al- 
lowed us to analyze weakly nonlinear AC response^, as 
well as the nonlinear DC conductance, order by order in 
the bias voltage. In contrast, so far the scattering matrix 
theory can be applied up to the second order nonlinearity 
and linear order AC. 

Using the concept of nonlinear characteristic potential 
developed in the response theoryEj, we have found that 
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the scattering matrix theory can actually be further de- 
veloped to apply to higher order nonlinear DC situations. 
In addition, recognizing that an AC transport problem 
requires the self-consistent solution of the Schrodinger 
equation coupled with Maxwell equations, we have found 
a way to derive both the external and internal responses 
in equal-footing within the scattering matrix approach. 
It is the purpose of this article to report these results. 
In particular, we shall start from the scattering matrix 
theory and formulate an approach which is appropriate 
for analyzing linear order dynamic conductance and the 
weakly nonlinear DC conductance beyond the second or- 
der. We emphasis the properties of electric current con- 
servation and gauge invariance, and these properties are 
maintained by considering electron-electron interactions. 
The approach developed here is particularly useful for 
nonlinear DC conductance calculation, and we shall ana- 
lyze the third order weakly nonlinear transport coefficient 
for a double-barrier tunneling diode. Since the approach 
presented here can be cast into a form which allows nu- 
merical computation, many further applications of it to 
complicated device structures can be envisioned. 

The rest of the article is organized as follows. In the 
next section we present the development of the formal- 
ism. Section III presents two applications of this formal- 
ism: the linear AC dynamic conductance and the third 
order nonlinear conductance. Finally a short conclusion 
is included in section IV. 



II. THEORETICAL FORMALISM 

In this section we briefly go through the formal devel- 
opment of our scattering matrix theory and concentrate 
more on the conceptually important physical quantities 
which will be needed. 

We start by writing the Hamiltonian of the system in 
the presence of an external time-dependent field as 

H = y^^(E am + eV a cosuit)al im (E am ,t)a am (E am ,t) 

am 

(i) 

where a) am is creation operator for a carrier in the incom- 
ing channel m in probe a, eV a cos cot is the shift of the 
electrochemical potential /i a away from the equilibrium 
state associated with fi eq , i.e., eV a cosu)t = fj, a — fi eq . 
The energy E am is a functional of the internal electrical 
potential landscape U(r,{V a }) which depends on V a in 
the low frequency regime. Potential U includes the in- 
ternal response to the external perturbation and it gen- 
erates such effect as the displacement current. In general 
U is also an explicit function of time (or of the AC fre- 
quency uj) as discussed in Ref. [li], but in this work we 
shall only be concerned with the dynamic conductance to 
first power of u> and for this case U is static. Note we have 
explicitly included U into the Hamiltonian which helps 



in dealing with both external and internal responses in 
equal-footing. The self-consistent nature of this Hamil- 
tonian is clear: U must be determined, in general, from 
the Maxwell equations where the charge density is ob- 
tained from solving the quantum mechanical problem of 
the Hamiltonian. 

Next let's consider the series expansion of the energy 
in terms of the potential landscape U, 
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where the operators O a is a spatial integration of the i- 
th order characteristic potential (see below) folded with 
the i-th order functional derivative of E am with respect 
to the potential landscape U(r). For instance the linear 
order operator, which is linear in voltage Vp, is given by 
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where 
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with uptr) = 9 gy* the linear order characteristic 
potentialo. The expressions for higher order operators 
are more difficult to write down in a general form, 
but they are proportional to the i-th power of the bias. In 
addition they can be easily determined after we formally 
obtain the transmission function and then applying the 
current conservation and gauge invariance to the results. 
Using Eq.(pl), the Hamiltonian now reads 



H = J2 
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The operators a am (E am , t) satisfies the equation of mo- 
tion 



in 
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which can be integrated because the time dependence of 
H is simple. For instance to linear order in voltage, we 
only need to use in the Hamiltonian and the result 
is 
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Its Fourier transform is given by 
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a a (E)= I dt& a (E a ,t)e iEt / n — a a (E) 

(6^) 2 [a a {E + 2hu)-2a a (E) 
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where we have suppressed the index m and a a is in a vec- 
tor form of the operators a am . In Eq.(^) the physics is 
transparent: a a (E ±?ilo) is just the one-photon sideband 
and a a (E ± 2Tioj) corresponds to the second harmonic 
generation. More tedious expressions can be obtained if 
higher order operators OW are included in the Hamilto- 
nian. 

To calculate the total electrical current, we shall apply 
the formula derived in Rcf. 15, which is exact up to linear 



order of uj and for larger frequency it is an approximation 
to a space-dependent expression of the current operator, 

I a (t) = J JdEdE' [at 

- hi (E)b a (E')\ x exp(*(£ - E')t/h) (7) 

where h a (E) is the operator which annihilates a carrier 
in the outgoing channel in probe a. The annihilation 
operator in the outgoing channel, b Q , is related to the 
annihilation operator in the incoming channel a a via the 
scattering matrix s Q ^: h a = ^2pS a pap where s a p is a 
function of energy E and a functional of the electric po- 
tential U(r, {V^}). Finally we comment that in evaluat- 
ing Eq. (0) we need to take a quantum statistical average 
of < at(E)ap(E') >= 5 a p6(E - E')f a (E) where f a (E) 
is the Fermi function of reservoir a. Because of the limi- 
tations of Eq. , our theory will be exact for transport 
coefficients linear in u> for AC situations. However, this is 
not a severe limitation for practical calculationsEj. 

One of the most important quantities of this theory 
is the determination of characteristic potential which ar- 
rives naturally. As discussed above, this quantity deter- 
mines the operators . Since the scattering matrix 
theory used here is exact to linear power of uj which is 
the order we shall work on, we only need to consider 
^-independent characteristic potentials. On the other 
hand, as we are interested in the weakly nonlinear coef- 
ficients, it is partial to consider higher order characteris- 
tic potentials^: ua 7 (r) = d 2 U(r)/dVsdV~ 1 , upsj( r ), etc.. 
For any physical quantity beyond the terms linear in uj or 
second order in voltage, including the second harmonic 
generation term (the term of ujV 2 ), these higher order 
characteristic potentials are necessary. 

We now discuss the solution of higher order character- 
istic potentials by explicitly carrying out the calculation 
of us-y- In the weakly nonlinear regime, the variation of 



the electric potential can be expanded in terms of the 
variation of the electrochemical potential dfi 



edU(r) = ^2 up{r)dfip + ^ ^ Upj(r)dfj.pd^ 
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where up is the characteristic potential, up-y (which is 
symmetric in (3 and 7) is the second order characteristic 
potential tensor, and (• • •) are higher order terms writ- 
ten in a similar fashion. Because we are only interested 
in AC transport to the first power of frequency w, the 
electrodynamics is solved by the Poisson equation, 



V 2 dU{v) = Aixe 2 dn(v) = 47re 2 V dn a (r) 



(9) 



where dn a is the variation of the charge density at con- 
tact a due to a change in electrochemical potential at 
that contact. There are two contributions to the charge 
density at contact a: the injected charge density due to 
the variation of the chemical potential at contact a, and 
the induced charge density dni n d, a due to the electro- 
static potential, hence 
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where dn a /dE is the injectivity which is the local den- 
sity of state at contact a and d 2 n a /dE 2 is the energy 
derivative of the injectivity. The induced charge den- 
sity involves Lindhard function. Using Thomas-Fermi 
approximation Eq.([Io|) takes a compact and simple form 
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From Eqs.(|), (§), and ([Ti"l), we obtain the equation sat- 
isfied by the second order characteristic potential tensor 
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Since all the quantities involved in this equation are 
known from the linear order calculation, up y can thus 
be determined. Similarly, order by order we can deter- 
mine higher order characteristic potentials from results 
obtained at lower orders, for instance the equation satis- 
fied by up-ys is found to be: 
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where the curl bracket {...} c stands for the cyclic permu- 
tation of indices (3, 7, and S. Note that if better models 
are needed to deal with the screening effect, the term with 
dn/dE on the left hand side of (H), @ and higher order 
equations is replaced by an integration over the appro- 
priate Luadhard function folded with the characteristic 
potentials. 

From Eq.(^j) we can derive several important sum 
rules on the characteristic potential tensor. If all the 
changes in the electrochemical potentials are the same, 
i.e., dfip = rf/i 7 = d/i, this corresponds to an overall 
shift of the electro-static potential edU = dfi. From 
this we have ^2pUp = 1, J2p-y u P~( = 0- Due to gauge 
invariance, Eq.(J8|) remains the same if dll, d/ip, and 
dfj,j are all shifted by the same amount. This leads to 
S/3 u 0i = S 7 u 0i = 0- Using Eq. (|l2|) we can confirm 
that these relations are indeed satisfied. Similar sum 
rules can be derived for higher order characteristic po- 
tentials. 

Let's summarize the scattering matrix theoretical pro- 
cedure. With the characteristic potential tensor calcu- 
lated, we explicitly derive the Hamiltonian in a series 
formEq. (|). The Hamiltonian determines the creation 
and annihilation operators via equation of motion Eq. 
(|J) and the scattering matrix s Q ^ . Finally, using Eq. (Q) 
we compute the electric current as a function of voltage. 

III. APPLICATIONS 

In the following we apply the scattering matrix for- 
malism developed in the last section to two examples: 
the linear order emittance and the third order DC non- 
linear conductance. The first example has been exam- 
ined by Biittiker and co-workersB, our result is in exact 
agreement with theirs. The second example has not been 
studied and we shall provide further numerical results for 
a resonant tunneling diode. 

A. Linear Dynamic Conductance 

The linear dynamic conductance (called emittance) is 
the transmission function of the terms proportional to Vp 
and ojVp in the electric current. From Eq. (Q) we expand 
every thing in these variables and obtain 

I a (u')=Y, 6 f I dE(-d E f){^App(a,E,E) 
p J 

e 2 + 

- —%U)'[& a pdES a p - {dES a p)& a p\} (14) 

where we have used notation A a p = App(a, E, E) and 
s a p = s a p(E, U). The transmission function A is defined 
in the usual form as 



App'(a,E,E',U) = l a 8apS a p> 

-sl {E,U) Sct p.{E',U) . (15) 

In deriving Eq. (|l~4|), we have used the fact that 
^2p App(a, E, E, U) = 0. All quantities such as A a p, 
s a p, and up — (dll (r) / dVp]e q are taken at equilibrium, 
i.e., at V a = 0. Using Eq.(Q), we separate the operator 

according to O^Oe — VpdE + J2-y V-ydv.,, thus ( |l~4|) can 
further be simplified to, 

I a (uj')= j dE(-d E f){J2^TA a p 
J p 

e 2 

~ ^h^Uj'^Vpls^pdESap ~ (dES a p)s a p] 
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From this result, we immediately realize that the first 
term on the right hand side is just the DC contribution 
to the electric current. From the second and third terms 
which are linear in lu' and Vp, we obtain the linear order 
emittance 

E a p = J dE(-d E f)-^—{[^ a pdES a p - {d E s a p)s a p} 

+ Yl^U dv f s ^ ~ ( 9 v>s a7 )s Q7 ]} . (17) 

7 

This result exactly agrees with the that obtained 
previouslyn. The first term of E a p describes the ex- 
ternal contribution to the AC current, while the second 
term is from internal response. They are obtained si- 
multaneously from the scattering matrix theory devek 
oped here. Finally, from the gauge invariance conditionEI 
ed E A a p + J2 7 d Vy A a p = 0, it is easy to show J^p E ap 

I ) . which is a direct consequence of O^ 1 ' in Eq. (|l4|). It is 
also easy to show that the electric current is conserved, 

Le -> J2 a E a0 = 0. 

B. Third Order DC Nonlinear Conductance 

The scattering matrix theory developed here can be 
applied to compute DC weakly nonlinear conductance to 
any order in bias. As an example we now calculate the 
third order DC nonlinear conductance G a p 7 s, which is 
defined by expanding the electric current in powers of 
voltage to the third power, 

I a = ^ G a pVp + G a p 1 VpV 1 

p p 1 

+ J2 G «0lS V V 7 V 6+ ■ ( 18 ) 

pjS 

Following the same procedure as above in deriving the 
linear emittance E a p, by expanding the electric current 
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Eq. (0) and other quantities to third order in bias, it is 
tedious but straightforward to derive 

e 3 f 

G a(ill) = — / dE(-d E f) [{dv„dv s A a(3 

+ edv y d E A af 35p S } c + e 2 d 2 E A a p5p 1 5 1 s\ . (19) 

Note that the second order characteristic potential tensor 
has been implicitly included in Eq. (|l9|) , because 



dv-,dv 5 A a/3 = d Vj 



5U(rt) dV s 
5 5A al 3 
5U(r 2 ) 8U( ri ) 
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(20) 



Again, we emphasis that other higher order nonlinear 
conductance can be calculated in a similar fashion. In 
general, the n-th order characteristic potential tensor 
is needed for the (n + l)th order nonlinear conduc- 
tance. Finally we point out that this result Eq.(|l9|) 
can in fact be obtained by expanding the following elec- 
tric current expression to the third order in voltage, 
4 = |EJ dEf(E -E F - eV )A a p(E, {V 7 }). 

In the following we calculate Gun from the general 
result of Eq. (|l]|) for a double barrier tunneling diode. 
For simplicity let's consider an one-dimensional double 
barrier tunneling system where the two barriers are 5- 
functions located at positions x — —a and x = a. The 
barrier strength is V\ and V<z respectively. When V\ = V2, 
this is a symmetric system, hence the second order non- 
linear conductance vanishes. In the symmetric case the 
first nonlinear coefficient comes from the third order, as 
specified by Eq. (|l9|). If we approximate the scattering 
matrix by the Brcit-Wigner forrrjj near a resonance en- 
ergy E r , s a p(E) ~ [6 a p - iy/TaTp/A], where T a is the 
decay width of barrier a, A = AE+iT/2 withT = Fi+r 2 
and AE = E — E r , we obtain a simple expression, 
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For the symmetric case, this expression reduces to 
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which is negative definite and has one minimum at AE — 
0. Because of the simple nature of the scattering matrix 
within the Breit-Wigner form- a general electric current 
expression has been obtainecE: I a = I a (Vp). We have 
thus calculated Gnu from this exact I-V relation, and 
it agrees exactly with the result (22) which comes from 
Eq. ©. 



Most practical transport problems can not be solved 
analytically. It is thus very important to be able to solve 
them numerically. Indeed, a distinct merit of the scat- 
tering matrix theory presented here is that it allows nu- 
merical computation, e.g. Eqs. (fT^), and ( [l9| ) can be nu- 
merically evaluated for explicit scattering potentials of a 
conductor. We only mention that the functional deriva- 
tives of the transmission function A a p with respect to 
the potential landscape U(r) as appeared in Eq. (|lS|), 
the potential derivatives, and the partial local density of 
states which is needed in the emittance calculation, can 
all be determined via the scattering Green's function. 

For the double-barrier diode just discussed, if we do 
not use the Breit-Wigner scattering form, Gnu can only 
be obtained numerically. For this systempthe Green's 
function G(x,x') can be calculated exactlyc3 thus from 
the Fisher-Lee relatioriiH we obtain the scattering matrix 
s a p(E) and hence the transmission function A a p from its 
definition (|l5|). To compute SAn/SU and the higher or- 
der functional derivative, we use Fisher-Lee relation and 
the factO that SG{xi,x 2 )/SU(x) = G(x 1 , x)G(x, x 2 ). 
Hence, 



fen 

6U(x) 



ihvG(x\,x)G(x, x\) 



(23) 



and 



+G(x 1 ,x)G(x,x')G(x',x 1 )} , (24) 
where v is the velocity of the particle. The energy and 



V 1 derivatives of Eq. <\lu) can be evaluated explici 



using the numerical procedures documented beforeQII 
Finally, the nonlinear characteristic potential u±i is ob- 
tained from Eq. (O). 
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FIG. 1. G1111 as a function of the scattering electron 
energy E for a double barrier tunneling diode with symmet- 
rical barriers. Solid line: numerical results by solving the full 
quantum scattering problem using Green's functions. Dashed 
curve: using the approximate Breit-Wigner form of the scat- 
tering matrix. The units of the quantities are set by % = 1, 
e = 1 and m = 1/2. 

The numerical result for Gnu as a function of the scat- 
tering electron energy E is plotted as the solid curve in 
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Fig. (1). Around the resonance energy E r w 2.0, Gnu 
takes large but negative values. Thus for electron Fermi 
energy in this range, the current-voltage characteristics 
will be nonlinear and this may result to negative differ- 
ential resistance. Such a behavior has important signifi- 

2 3. Notice that the two dips 



cance for practical purpose: 
of Gnu are not exactly the same, such an asymmet- 
rical behavior has been observed in other quantities!!^. 
For comparison we also plotted (dashed line) the result 
from the Breit-Wigner approximation of the scattering 
matrix, Eq. (p2|). While the negative nature of Gnu 
and the overall magnitude are similar to the numeri- 
cal result, the Breit-Wigner result shows only one dip. 
This inconsistency is completely due to the simple form 
of the Breit-Wigner approximation as it gives a space- 
independent and constant characteristic potential. The 
accurate solution using the Green's function generates 
space-dependence of various quantities. 



IV. SUMMARY 

In summary, we have extended the scattering matrix 
theory which is now appropriate to analyzing linear dy- 
namic conductance to first order in frequency, and weakly 
nonlinear DC conductance order by order in external 
bias. The crucial ingredient of this development is the 
characteristic potential at weakly nonlinear orders and 
these potentials appear naturally from the self-consistent 
Hamiltonian. The theory is current conserving and gauge 
invariant. The physical quantities involved in this theory 
are numerically calculable, hence the present approach 
can be used to conductors with complicated scattering 
potential landscape for quantitative predictions. The 
formal connection of the scattering matrix theory and 
the response theory, at the weakly nonlinear AC level, 
has been clarified in our recpat workE-i While the re- 
sponse theory is very generalliil and can be used to ana- 
lyze weakly nonlinear DC and AC transport order by or- 
der in bias, and for AC case to all orders of frequency, we 
expect that the scattering matrix theory should be able 
to do the same. This paper partially fulfills this expecta- 
tion by extending the scattering matrix theory to higher 
orders of nonlinearity. We point out that to go beyond 
the linear frequency, the expression for electric current, 
Eq. (R), should be extended. In addition we should use 
the Helmhotz equation (in Lorentz gauge) for the electro- 
dynamics instead of the Poisson equation. Finally, it is 
important to note that the theoretical approach of scat- 
tering matrix developed here is appropriate to transport 
problems near equilibrium. Far from equilibrium, one 
may employ the Keldysh Green's functionsEd. 
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